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The quantum thermodynamic behavior of small systems is investigated 
in presence of finite quantum dissipation. We consider the archetype cases 
of a damped harmonic oscillator and a free quantum Brownian particle. A 
main finding is that quantum dissipation helps to ensure the validity of the 
Third Law. For the quantum oscillator, finite damping replaces the zero- 
coupling result of an exponential suppression of the specific heat at low 
temperatures by a power-law behavior. Rather intriguing is the behavior 
of the free quantum Brownian particle. In this case, quantum dissipation 
is able to restore the Third Law: Instead of being constant down to zero 
temperature, the specific heat now vanishes proportional to temperature 
with an amplitude that is inversely proportional to the ohmic dissipation 
strength. A distinct subtlety of finite quantum dissipation is the result 
that the various thermodynamic functions of the sub-system do not only 
depend on the dissipation strength but depend as well on the prescription 
employed in their definition. 

PACS numbers: 05.70.-a, 05.30.-d, 05.40.-a, 05.40.Jc 



1. Introduction 

The development of the theory of Brownian motion played a pivotal 
role - and continues to do so - in the development of statistical mechanics 
and thermodynamics [1,2]. Thermodynamics together with relativity and 
quantum theory form three pillars on which much of the entire structure 
of physics rests. Tampering with the axioms in either of those theories is 
not a good idea; doing so may well lead to contradictions with the other 
theories. In particular, the field of thermodynamics bears consequences for 
many branches of physics. Its four laws are well-known [3,4]: the zeroth law 



* This work is dedicated to our colleague and friend, Professor Peter Talkner, on the 
occasion of his 60-th birthday. 



CD 



2 



Peter Hanggi and Gert-Ludwig Ingold 



guarantees that states of thermal equilibrium exist which can be character- 
ized by a temperature T. The first law provides a balance among the various 
contributions that make up the internal energy of a system while the second 
law introduces the concept of thermodynamic entropy S, which notably is 
extensive and never decreases for a closed physical system. In addition, the 
second law tells us that there exists an absolute zero of temperature. 

The Third Law is attributed to Walther Hermann Nernst (1864-1941) 
and arose as the result of his seminal idea - being guided by his critical 
analysis of chemical and electrochemical reactions at lower temperatures - 
that at low temperatures there occurs for isothermal processes a perfect 
correspondence between the enthalpy and the Gibbs free energy. Thereby, 
the approximate rule hypothesized earlier by Marcelin Berthelot and Julius 
Thomson, becomes a "law" at zero temperature. Nernst announced this 
result already in his lectures in 1905, terming it "mein Warmesatz" (my 
law of heat) [5,6]. He took this result even further: He also studied how 
fast the difference between the changes in the enthalpy AH and the Gibbs 
free energy AG, i.e. AH — AG tends to zero [7]. In fact, this difference 
vanishes faster than linear in temperature implying that the change of en- 
tropy itself must vanish at absolute zero. This in turn implies identical, 
generally vanishing initial slopes for the corresponding quantities AH and 
AG as a function of temperature. An elucidating account of the history of 
the Third Law and the controversies surrounding its acceptance is presented 
in the books by Dugdale [8] and by Wilks [9]. In its strict form, as given by 
Max Planck [10], the Third Law reads: The entropy s = S/N per particle 
approaches at absolute zero a constant value so that possibly depends on 
the chemical composition of the system. This limiting entropy constant so 
can generally be set equal to zero. 

The Third Law carries prominent consequences for quantum mechanics 
and the field of low-temperature physics. First, the fact that at absolute 
zero temperature the isotherm coincides with the isentrope (adiabat) imme- 
diately implies that this absolute zero temperature is unattainable by use 
of a sequence of isothermal and adiabatic reversible operations [11]. There- 
fore, it has the consequence that the efficiency of a Carnot engine, i.e. a heat 
engine that cyclically operates between two heat baths of different tempera- 
tures which are never brought into contact with each other, can never reach 
100% for any finite upper temperature. Moreover, the constant value of 
the entropy at absolute zero is given by the degeneracy g of the iV-particle 
system in the corresponding quantum ground state, i.e. S(T = 0) = &b In g, 
where k-g is the Boltzmann constant. The limiting value of the intensive 
quantity sq = S(T = 0)/N in the thermodynamic limit of particle number 
N — > oo will typically be zero, so long as the degeneracy g = g(N) does not 
grow with N faster than exponentially [12]. A well-known exception is the 
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case of noninteracting, independent particles carrying a non- vanishing spin 
/, yielding sq = k-Q ln(2I+l) for the limiting entropy per particle. Moreover, 
the Third Law also implies that thermal quantities such as specific heats, 
the isobar thermal coefficient of expansion, the isochor coefficient of tension, 
etc., all approach zero as T — > 0. Likewise, the (magnetic) susceptibility 
becomes constant as T — > 0, so that the classic Curie law must loose its 
validity at very small temperatures. 

Are there known exceptions of systems not obeying the Third Law? 
It is known that many classical systems do not obey the Third Law. In 
particular, noninteracting classical systems with their constant values for 
the specific heat clearly violate the Third Law. A well-known case is 
the classical ideal gas for which the entropy S assumes the form S = 
N[cy InT + fee ln(V/N) + a], where V is the volume and a denotes the en- 
tropy constant. It clearly does not fulfill the Third Law because it diverges 
logarithmically with temperature T for a constant specific heat cy. Even 
when we use for cy the physically relevant low-temperature result, namely 
that quantum mechanically the specific heat cy vanishes faster than InT, 
we still find a dependence on the particle density which is not compatible 
with the formulation of Planck. This observation that the classical gas does 
not obey the Third Law led Nernst to speculate that the classical gases must 
undergo a "degeneracy", which has been resolved with the quantum statis- 
tics for the ideal Fermi gas and the ideal Bose-Einstein gas, which indeed 
do obey the Third Law in the strict formulation by Planck. 

Are there yet other remaining open problems with the Third Law? Ap- 
parent difficulties with the Third Law occur for metastable states that do 
not necessarily guarantee the sufficiently fast relaxation within a finite time 
scale towards thermal equilibrium, the latter being a prerequisite for the va- 
lidity of the Third Law. In this context, glasses provide a system class that 
can provide detectable deviations from the Third Law at low temperatures 
which likely are the result of frozen-in ordered excited states that have not 
yet fully relaxed. According to common wisdom the known deviations from 
the Third Law will all be cured by quantum mechanics, quantum statistics, 
and interactions among particles. 

In the following we shall investigate the thermodynamic low tempera- 
ture properties for open quantum systems that are coupled to a heat bath of 
finite dissipation strength. Because both, the thermodynamics of classical 
open systems and the quantum statistical mechanics of open systems are 
strictly valid only for systems that are only infinitely weakly coupled to a 
bath, it is a priori not obvious how the Third Law lives up to the sub-system 
in presence of finite quantum dissipation [13-22]. The effect of the finite 
coupling of a sub-system to an environment in fact induces several subtleties 
for quantum Brownian motion [22]. For example, the equilibrium density 
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matrix is no longer given by its standard canonical form p$ = exp{—(3H^)/Z 
where H$ denotes the system Hamiltonian and Z = Tr[exp(— /Wis)] is the 
partition function. Therefore, in clear contrast to the classical case with 
dissipation, this reduced density matrix becomes also a function of the in- 
teraction strength with the environment. This being so, taking a gas of 
free independent quantum Brownian particles that are coupled to a heat 
bath with finite dissipation strength, one may speculate that the role of the 
interactions of the free particle with the abundant bath degrees of freedom 
will be sufficient to cure the shortcomings stemming from a classical gas of 
free Brownian particles. 

2. Quantum harmonic oscillator revisited 

To start out, we first recall the standard results for a single harmonic 
quantum oscillator that is infinitely weakly coupled to a bath that estab- 
lishes the temperature T. This situation is reminiscent of the famous treat- 
ment of the specific heat of a solid by Albert Einstein in 1907 [23], where he 
found an exponential suppression of the specific heat as T — > 0. This find- 
ing impressed the Berlin school so immensely, in particular Nernst and his 
collaborators (who in 1910 experimentally confirmed this salient first pre- 
diction of quantum theory) , that Nernst together with Planck were able to 
bring the "new Copernicus" [5,6] into the exclusive circle of Berlin physicists 
in 1913. 

2. 1 . Partition function and entropy 

Let us consider an oscillator degree of freedom of mass M and force 
constant /, i.e. its Hamiltonian H$ reads 

H s = m + \^- <»> 

In terms of the angular frequency ujq = f /M the quantum mechanical en- 
ergy eigenvalues read E n = (n + ^)hu>o, yielding for the partition function 
Z the well-known expression 

Z = V e~P En = - (2) 

^ 2sinh[?i/3wo/2] K ' 

where (3 = l/k^T is the inverse temperature. 
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Using familiar relations we find that the internal energy E reads 

E=-— \n(Z) (3) 



+ 



flUJn 



2 ' ex.p(h(3u!o) 
and, correspondingly, the entropy is given by 
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In (l — exp(— h(3tdo)) 



exp(/i/?u;o) — 1 

For low temperatures, hfluiQ S> 1, the entropy approaches zero like 



5 



exp 



Hun 



T L \ k B T / 
The specific heat can now be derived either from as 



BE 



or from © as 



C = T d J: = -P? S 



dT 'dp' 
In both cases, one obtains for the specific heat 

2smh(h/3io /2) J 
Its low-temperature behavior 

C = k B 



C = k B 



[m exp 



flUJn 



k B T 



(4) 

(5) 
(6) 

(7) 

(8) 
(9) 

(10) 
(11) 



is not analytic in temperature and corresponds to Einstein's result for the 
low-temperature behavior of the specific heat of a solid [23]. For high tem- 
peratures one finds 

\ _ 1 (^1 2 



C = k B 



12 \k B T 



+ 0{T- 



(12) 



As this result shows, the specific heat for a free particle cannot simply be 
obtained by taking the limit uq — > of the harmonic oscillator. Such a pro- 
cedure will not properly account for the reduced number of degrees of free- 
dom which within the equipartition theorem will lead to a high-temperature 
specific heat of only C = k B /2 for the free particle. 
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3. Quantum harmonic oscillator: The role of quantum dissipation 

3.1. Harmonic oscillator coupled to an environment 

We now couple the harmonic oscillator of the previous section to an en- 
vironment consisting of an infinite number of harmonic oscillators forming a 
heat bath. In contrast to the previous section, the coupling strength will not 
be kept negligible here. In addition, system and bath are infinitely weakly 
coupled to a superbath which has the purpose to provide the temperature 
T. 

The total Hamiltonian H does not need to account for the superbath 
and therefore consists of three parts [15, 16, 18, 22, 24, 25] 

H = H S + H B + H SB (13) 

where H§ is given by 1(1)). the bath Hamiltonian reads 

+ (14) 

and the coupling is bilinear in the coordinates 

°° c 2 

tfsB = -4£ C ^ + S 2 £^2- (15) 
i=l 1 

We note that quantum systems that are coupled to an environment of 
finite strength are rarely exactly solvable. The dissipative quantum oscil- 
lator becomes exactly solvable with its bilinear coupling to a bath because 
of the inherent quadratic structure of the total Hamiltonian in ()13|) . This 
fact holds true even for the case of time-dependent, parametrically driven 
dissipative quantum harmonic systems [26,27]. 

In order to describe the influence of the environment on the system 
oscillator, it is sufficient to know the spectral density of bath oscillators 
defined by [13, 15, 16, 19] 

oo 2 

J( w ) = 7r£— i-*(a,-u, f ). (16) 
i=i 

For later purposes, we introduce the Laplace transform of the damping 
kernel, which generally depends on frequency, thereby causing memory- 
friction [28,29], i.e., 

If 00 dwJ(w) 2z 

The important special case of strictly ohmic dissipation is characterized 
by J(cj) = M^fu and 7(2;) = 7 which leads to a memoryless damping of 
strength 7. 
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3.2. Specific heat of a damped harmonic oscillator 

We next discuss the specific heat of the damped harmonic oscillator by 
following two routes. First, we start from the energy E and employ the 
common relation in ([%]). As an alternative route we shall in Section 13.31 
determine the entropy from the partition function by means of ® from 
which allows one to evaluate the specific heat in the case of strictly 
ohmic damping. 

The energy of the damped harmonic oscillator is given by 

<*H§ + ^oV> (is) 

where the expectation value of an operator 0$ acting in the Hilbert space of 
the system is defined with respect to the canonical density matrix of system 
plus environment as 

For ohmic damping, the second moments of position and momentum can 
be expressed as [15,30,31] 

(q 2 ) = ^fo(T) (20) 



and 

(p 2 )=Mhf 2 (T) (21) 
where we have introduced a temperature-dependent function 

2tt (lo 2 - to 2 ) 2 + 7 2 w 



For n = 2, i.e. when evaluating (p 2 ), the integrand decreases only with 1/lo 
and a finite value can only be obtained by introducing a high-frequency cut- 
off in the damping kernel ■j(z). However, this divergent term gives rise only 
to a temperature-independent contribution to {p 2 ) and thus to the energy 
(|18|). When evaluating the specific heat according to |HJ), this constant term 
will disappear and a finite result is obtained even for ohmic damping. After 
some algebra, one finds for the specific heat 

£ = 1 - ^ + W (1 + A+) + A 2 1/ (1 + A_) (23) 
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where 



2ir 



7 
2uj 




(24) 



and vp'(z) is the trigamma function. At low temperatures, the specific heat 
thus assumes the form 



C_ 

k B 



7r 7 k B T 4-7T 3 7 



3 ujq Tiujq 15 loq 



- 



(25) 



This result differs significantly from the expression qiljl in the absence of 
dissipation. While in the latter case, the presence of an energy gap led to an 
exponential suppression of the specific heat, we now find a linear increase 
with temperature. This behavior indicates the existence of a finite density 
of states even at small excitation energies [32]. Even at high temperatures 
the effect of dissipation can be detected, albeit in a less spectacular manner. 
The leading correction in the high-temperature expansion reads 



C_ 

k B 



h'j ^ 2 (7 



+ 



2wg) 



+ 0(T~ 



(26) 



2irk B T ' 24(fc B T) 2 

Thus, the leading correction depends on the damping strength 7. This 
finding is in clear contrast to the behavior of the quantum escape rate 
[16,29]: There, the leading quantum correction to the escape rate always 
enhances the classical result and is independent of the dissipation strength. 



3.3. Energy versus partition function for a damped harmonic oscillator 

Another prescription to obtain the specific heat starts out from the 
canonical partition function 

Tr [exp(-/3H)] 



Z 



Tt b [exp(-pH B ) 



(27) 



where Tr^ denotes the partial trace in the Hilbert space of the bath. In the 
absence of a system-bath coupling, this expression would correspond to the 
partition function of the system alone. For a damped harmonic oscillator, 
the partition function becomes [33] 



1 



n 



(28) 



with the Matsubara frequencies v n = 2irn/hf3. In view of the divergence 
for strictly ohmic damping mentioned above, we allow here for a possible 
frequency dependence of 7. 
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Following the standard procedure of statistical mechanics, we can obtain 



the energy from the partition function by means of 



<£>z=-|mz). 



Inserting (|28|). one obtains 



(Eh 



1 



71=1 



2lvI + Vn^jVn) - vjl'^n) 



which in general differs from the expression 



(E) 



1 

1 



2ug + t/ n 7(i/ TO ) 



(29) 



(30) 



(31) 



obtained by evaluating the integral (|22|) by residues. The only exception is 
the special case of strictly ohmic damping where j(v n ) = 7 is constant. 

This generally non-vanishing difference does not come as a surprise if 
one only takes a closer look at the partition function (|27|): The evaluation 
of (JH yields 



(E) z = (H) - (H b )b 

= (E) + [{Hsb) + (Hb) 



(Hb)b] 



(32) 



where the index "B" denotes an average with respect to the bath Hamilto- 
nian only. This result differs from the energy (E) = (H$) by the term 
in the brackets which, generally, vanishes only in the absence of a system- 
bath coupling. The coincidence between (E) and (E)z for the harmonic 
oscillator subject to strictly ohmic damping should therefore be considered 
as exceptional. 

Nevertheless, we briefly sketch how one would obtain the specific heat 
from the partition function because this will give us as a by-product an 
expression for the entropy of the damped harmonic oscillator. For strictly 
ohmic damping the product 1)28(1 does not converge and in principle a high- 
frequency cutoff for 7 should be introduced. However, this divergence can 
again be traced back to an infinite energy shift due to the environmental 
coupling. We may shift the energy by an arbitrary amount A by multiplying 
the partition function by exp(— (3A) without changing the entropy or the 
specific heat. After performing an appropriate energy shift, we arrive at an 
expression of the partition function valid even for strictly ohmic damping 
[34], reading 



A ( 2vr 

h(3u>Q \h(3u>o 



r(i + A+)r(i + A_) 



(33) 
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where T(z) is the gamma function, A± are defined in (|24j). and A is a 
constant whose precise value is irrelevant for the following. By virtue of JSJ 
we obtain for the entropy the result 



S = k B 



1 - \n(h/3u ) + 



nth 

2tt 



(34) 



where we introduced the abbreviation 



g(z) = ln[r(l + *)] - zij){l + z) (35) 

with the digamma function ip{z) = T'(z)/T(z). In the absence of damping, 
i.e. 7 = 0, this reproduces the result © for the entropy of an uncoupled 
harmonic oscillator. 

For very low temperatures, the entropy (|34|) vanishes like 

s= * :L k^T +0{T3) 
6 ujq nujQ 

as required by the Third Law of thermodynamics. By means of ©, the 
expression (|23|) for the specific heat is recovered identically. 

The low-temperature behavior (|36j) is in agreement with the expression 
derived by Ford and O'Connell on the basis of the free energy [35]. For the 
damped harmonic oscillator these authors found that the entropy vanishes 
in the limit of zero temperature also for more general forms of the bath 
density of states. 



4. Free quantum Brownian motion coupled to a heat bath: Is the 
Third Law obeyed at zero temperature? 

According to the equipartition theorem the specific heat of a free parti- 
cle is /cb/2. In the limit of an infinite "box" this represents even the correct 
quantum value, because the classical and the quantum partition function 
become equal. In view of the fact that the specific heat of an ideal gas thus 
remains non-zero down to the lowest temperatures, one becomes curious to 
investigate the specific heat of a free particle coupled to an environment 
when the coupling strength is not assumed to vanish: Does quantum dissi- 
pation help to restore the Third Law? 

Free quantum Brownian motion has been addressed in earlier work 
[36-39] wherein the main focus centered on the role of free quantum diffu- 
sion [15]. Interestingly enough, for ohmic dissipation the quantum diffusion 
remains "classical" , being proportional to time t, except at zero temperature 
itself, where one finds a logarithmic behavior in time t [36-39]. At finite 
temperatures this quantum behavior is observable at intermediate times 
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only [36,40]. One is thus tempted to conclude that finite quantum dissipa- 
tion will not be sufficient to cure the classical behavior for the specific heat. 
Therefore, we shall next investigate the behavior of the specific heat for free 
quantum Brownian motion in closer detail. 

As already mentioned at the end of Section 12.11 simply taking the limit 
ujq — > in the results for the damped harmonic oscillator is not without 
problems. We therefore start our calculation from the energy 



W 2M 2/3 



(37) 



which can be obtained from the second moment of momentum {p 2 } by evalu- 
ating the integral (|22|) by residues. Proceeding as in Section l3~2l one derives 
the specific heat of the free particle in the presence of strictly ohmic damping 

k B 2 2vr V 27r / V 27T J 

which for either T — > oo or 7 — > yields C = /cb/2, as expected from the 
equipartition theorem. On the other hand, for low temperatures the specific 
heat tends to zero as 

kb 3 nj 15 \ nj J 

in agreement with the Third Law of thermodynamics. The specific heat (|38|) 
together with its linear low-temperature behavior are depicted in the main 
part of Fig. ^ as full and dashed line, respectively. Temperatures krT 3> ^7 
much larger than the damping strength are required in order to restore the 
classical result C/k-Q = 1/2. 

As the damping strength already serves to set the temperature scale, it 
is instructive to introduce a cutoff in the density of states (|16|) of the bath 
oscillators in order to study how a reduction of the environmental influence 
changes the dependence of the specific heat on temperature. To this end, we 
introduce a high-frequency cutoff ujd by choosing the Drude model where 

l{z) = ■ (40) 

This environment leads to a specific heat 
C 1 n/3 7 1 



[z + ifj'(l + z + )-z^/(l + z^)] , (41) 
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Fig. 1. The specific heat C of a free quantum Brownian particle is shown as a 
function of the temperature T for strictly ohmic friction of strength 7. The dashed 
line indicates the linear low-temperature behavior. In the inset, the modification 
due to a finite cutoff frequency wd is depicted for different cutoff scales, ivu/l — 
0.01, 0.1, 1, and 00 from the upper to the lower full line. The dashed line indicates 
again the linear low-temperature behavior. 

where we have introduced the abbreviations 



Independently of the value of the cutoff frequency wd, the specific heat 
will go to zero linearly as stated in 1)39(1 . However, for temperatures larger 
than the cutoff frequency, the suppression of the specific heat due to the 
environmental coupling will be ineffective. A reduction of the environment 
thus tends to restore the classical value of fee/2 for the specific heat. This 
effect is depicted in the inset in Fig. ^ where the cutoff frequency takes 
the values = 0.01,0.1,1 and 00 from the upper to the lower curve. 

The dashed curve represents the linear low-temperature behavior which still 
dominates at temperatures below Tiuj^/kB- In contrast to many phenomena 
in quantum dissipation, here an increase of the coupling to the environment 
does not render the system more classical. On the contrary, a stronger 
environmental coupling makes the dissipative quantum system behave more 
quantum mechanically and thus helps to ensure the validity of the Third Law 
of thermodynamics, i.e. the vanishing of the specific heat with decreasing 
temperatures. 




(42) 
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5. Conclusions 

With this work we have explored the behavior of the specific heat for a 
quantum system that is coupled to a heat bath with finite coupling strength. 
Our findings are contrasted with the Third Law of thermodynamics which 
generically predicts a vanishing of the specific heat at low temperatures. For 
a harmonic oscillator the presence of quantum dissipation changes the well- 
known Einstein-like behavior of an exponentially fast approach towards zero 
specific heat into a power-law behavior with a slope that increases with in- 
creasing coupling strength. Even more intriguing is the behavior for a freely 
moving quantum particle: While the quantum treatment in absence of dis- 
sipation simply coincides with the classical behavior, i.e. the specific heat 
takes a constant value C = ks/2, the role of finite quantum dissipation is 
able to restore the Third Law, yielding a leading linear temperature depen- 
dence. Quite counterintuitively, its approach to the classical value occurs 
the faster the weaker is the dissipation strength. 

In contrast to common quantum statistical mechanics which intrinsically 
is based on a vanishingly small coupling to the environment, the finite cou- 
pling strength between the sub-system and the bath causes some subtleties 
that must be recognized. As made explicit in Section 13.31 the thermody- 
namic quantities depend on the procedure invoked in their definition: the 
commonly used expression based on the partition function provides results 
that generally do not agree with the result obtained from the correspond- 
ing quantum expectation value. Interestingly enough, in the strict ohmic 
limit (i.e. in the absence of a high-frequency cutoff) the specific heat for the 
damped harmonic oscillator does not depend on the prescription employed. 
For the case of memory friction, however, where a finite cutoff frequency 
is present, the thermodynamic quantities depend both on the value of this 
cutoff and the prescription used in their evaluation. 

The results obtained for the low-temperature behavior of the specific 
heat of simple quantum systems are not only of academic interest, but 
may turn out to be relevant for experiments in nanoscience where one tests 
the quantum thermodynamics of small systems [41] that are coupled to an 
environment with a finite coupling strength. 
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